---------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/log/Table_A3_Fig5_Results.log
  log type:  text
 opened on:  15 May 2024, 10:56:10

. eststo clear

. use "$rootpath/processed_data/Estimation_Data.dta", clear

. drop if key_estimation_sample==0
(273,947 observations deleted)

. label variable lead_5 "-4"

. label variable lead_4 "-3"

. label variable lead_3 "-2"

. label variable lead_2 "-1"

. label variable lead_1 "0"

. label variable lag_1 "1"

. label variable lag_2 "2"

. label variable lag_3 "3"

. label variable lag_4 "4"

. label variable lag_5 "5"

. sort facility_id year

. egen num_obs=nvals(year), by(facility_id)

. 
. rename lead_1 lead_0

. //facility-level errors
. reghdfe  log_damage55 lead_5 lead_4 lead_3 lead_2 lead_0 lag_* if num_obs==10&min_year==2008, absorb(facility_id naics#year) vce(cluster facility_id)
(dropped 1414 singleton observations)
note: lead_0 is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
(MWFE estimator converged in 205 iterations)
note: lead_0 omitted because of collinearity

HDFE Linear regression                            Number of obs   =    313,926
Absorbing 2 HDFE groups                           F(   9,  31475) =       7.49
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9444
                                                  Adj R-squared   =     0.9363
                                                  Within R-sq.    =     0.0005
Number of clusters (facility_id) =     31,476     Root MSE        =     0.4192

                       (Std. Err. adjusted for 31,476 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      lead_5 |  -.0101698    .030879    -0.33   0.742    -.0706939    .0503543
      lead_4 |   .0161482   .0273878     0.59   0.555    -.0375329    .0698293
      lead_3 |   .0251618   .0273605     0.92   0.358    -.0284658    .0787894
      lead_2 |   .0115403    .023941     0.48   0.630     -.035385    .0584656
      lead_0 |          0  (omitted)
       lag_1 |  -.0010444   .0141349    -0.07   0.941    -.0287494    .0266606
       lag_2 |  -.0338259   .0236071    -1.43   0.152    -.0800967    .0124449
       lag_3 |  -.1611297    .027961    -5.76   0.000    -.2159344   -.1063251
       lag_4 |  -.0828922   .0300817    -2.76   0.006    -.1418535   -.0239309
       lag_5 |   -.123768   .0308232    -4.02   0.000    -.1841827   -.0633532
       _cons |   1.321442   .0060274   219.24   0.000     1.309628    1.333255
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     31476       31476           0    *|
   naics#year |      8154           1        8153     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store event_study44

. coefplot event_study44, level(99) vert keep( lead_5 lead_4 lead_3 lead_2 lead_0 lag_1 lag_2 lag_3 lag_4 lag_5) yline(0) nokey omitted ytitle("Impact on
>  Inverse Hyperbolic Sine (Damages)") xtitle("Years from Policy Change") ciopts(recast(rcap)) graphregion(color(white))

. graph export "$rootpath/table_graph/Event_Study_Damages.pdf", replace
(file /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/Event_Study_Damages.pdf written in PDF format)

. 
. //state-level errors
. reghdfe  log_damage55 lead_5 lead_4 lead_3 lead_2 lead_0 lag_* if num_obs==10&min_year==2008, absorb(facility_id naics#year) vce(cluster state_number)
(dropped 1414 singleton observations)
note: lead_0 is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
(MWFE estimator converged in 205 iterations)
note: lead_0 omitted because of collinearity

HDFE Linear regression                            Number of obs   =    313,926
Absorbing 2 HDFE groups                           F(   9,     47) =       3.59
Statistics robust to heteroskedasticity           Prob > F        =     0.0018
                                                  R-squared       =     0.9444
                                                  Adj R-squared   =     0.9363
                                                  Within R-sq.    =     0.0005
Number of clusters (state_number) =         48    Root MSE        =     0.4192

                          (Std. Err. adjusted for 48 clusters in state_number)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      lead_5 |  -.0101698   .0466907    -0.22   0.829    -.1040994    .0837597
      lead_4 |   .0161482   .0560346     0.29   0.774    -.0965789    .1288753
      lead_3 |   .0251618   .0477104     0.53   0.600    -.0708192    .1211427
      lead_2 |   .0115403   .0435027     0.27   0.792    -.0759758    .0990564
      lead_0 |          0  (omitted)
       lag_1 |  -.0010444    .036757    -0.03   0.977    -.0749899    .0729011
       lag_2 |  -.0338259     .06103    -0.55   0.582    -.1566025    .0889507
       lag_3 |  -.1611297   .0570047    -2.83   0.007    -.2758084    -.046451
       lag_4 |  -.0828922   .0647682    -1.28   0.207     -.213189    .0474046
       lag_5 |   -.123768   .0670302    -1.85   0.071    -.2586154    .0110794
       _cons |   1.321442   .0108078   122.27   0.000     1.299699    1.343184
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     31476       31476           0    *|
   naics#year |      8154           1        8153     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. 
. //facility-level se
. reghdfe  log_emissions4 lead_5 lead_4 lead_3 lead_2 lead_0 lag_* if num_obs==10&min_year==2008, absorb(facility_id naics#year) vce(cluster facility_id)
(dropped 1414 singleton observations)
note: lead_0 is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
(MWFE estimator converged in 203 iterations)
note: lead_0 omitted because of collinearity

HDFE Linear regression                            Number of obs   =    313,926
Absorbing 2 HDFE groups                           F(   9,  31475) =      14.43
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9473
                                                  Adj R-squared   =     0.9397
                                                  Within R-sq.    =     0.0009
Number of clusters (facility_id) =     31,476     Root MSE        =     0.3493

                       (Std. Err. adjusted for 31,476 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      lead_5 |   .0410573   .0264142     1.55   0.120    -.0107155    .0928301
      lead_4 |   .0571463   .0232722     2.46   0.014      .011532    .1027607
      lead_3 |   .0542545   .0232088     2.34   0.019     .0087643    .0997447
      lead_2 |    .000727   .0204616     0.04   0.972    -.0393786    .0408326
      lead_0 |          0  (omitted)
       lag_1 |  -.0193149   .0122987    -1.57   0.116    -.0434208     .004791
       lag_2 |   .0150585   .0229123     0.66   0.511    -.0298504    .0599675
       lag_3 |  -.1690039   .0241435    -7.00   0.000    -.2163262   -.1216817
       lag_4 |  -.0142423    .027678    -0.51   0.607    -.0684924    .0400077
       lag_5 |  -.1216359   .0262117    -4.64   0.000    -.1730119     -.07026
       _cons |   1.046349   .0051786   202.05   0.000     1.036199    1.056499
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     31476       31476           0    *|
   naics#year |      8154           1        8153     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store event_studyem4

. coefplot event_studyem4, level(99) vert keep( lead_5 lead_4 lead_3 lead_2 lead_0 lag_1 lag_2 lag_3 lag_4 lag_5) yline(0) nokey omitted ytitle("Impact o
> n Inverse Hyperbolic Sine (Emissions)") xtitle("Years from Policy Change") ciopts(recast(rcap)) graphregion(color(white))

. graph export "$rootpath/table_graph/Event_Study_Emissions.pdf", replace
(file /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/Event_Study_Emissions.pdf written in PDF forma
> t)

. 
. //state-level se
. reghdfe  log_emissions4 lead_5 lead_4 lead_3 lead_2 lead_0 lag_* if num_obs==10&min_year==2008, absorb(facility_id naics#year) vce(cluster state_number
> )
(dropped 1414 singleton observations)
note: lead_0 is probably collinear with the fixed effects (all partialled-out values are close to zero; tol = 1.0e-09)
(MWFE estimator converged in 203 iterations)
note: lead_0 omitted because of collinearity

HDFE Linear regression                            Number of obs   =    313,926
Absorbing 2 HDFE groups                           F(   9,     47) =       4.56
Statistics robust to heteroskedasticity           Prob > F        =     0.0002
                                                  R-squared       =     0.9473
                                                  Adj R-squared   =     0.9397
                                                  Within R-sq.    =     0.0009
Number of clusters (state_number) =         48    Root MSE        =     0.3493

                          (Std. Err. adjusted for 48 clusters in state_number)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      lead_5 |   .0410573    .039212     1.05   0.300     -.037827    .1199416
      lead_4 |   .0571463   .0378855     1.51   0.138    -.0190695    .1333621
      lead_3 |   .0542545   .0387892     1.40   0.168    -.0237793    .1322884
      lead_2 |    .000727   .0441723     0.02   0.987    -.0881361    .0895901
      lead_0 |          0  (omitted)
       lag_1 |  -.0193149   .0217648    -0.89   0.379       -.0631    .0244702
       lag_2 |   .0150585   .0779194     0.19   0.848     -.141695    .1718121
       lag_3 |  -.1690039   .0556491    -3.04   0.004    -.2809555   -.0570524
       lag_4 |  -.0142423   .0795318    -0.18   0.859    -.1742397     .145755
       lag_5 |  -.1216359   .0545713    -2.23   0.031    -.2314192   -.0118527
       _cons |   1.046349   .0076255   137.22   0.000     1.031008     1.06169
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     31476       31476           0    *|
   naics#year |      8154           1        8153     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. 
. esttab event_study44 event_studyem4 using "$rootpath/table_graph/EVENT_STUDY_Results_A3.tex",  star( * 0.1 ** 0.05 *** 0.01) cells(b(star fmt(3)) se(pa
> r fmt(3))) keep(lead_5 lead_4 lead_3 lead_2 lag_*) scalars(N F) ar2 replace 
(note: file /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/EVENT_STUDY_Results_A3.tex not found)
(output written to /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/EVENT_STUDY_Results_A3.tex)

. log close
      name:  <unnamed>
       log:  /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/log/Table_A3_Fig5_Results.log
  log type:  text
 closed on:  15 May 2024, 10:58:53
---------------------------------------------------------------------------------------------------------------------------------------------------------
